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Abstract. Dose calculation for radiotherapy with protons and heavier ions deals 
with a large volume of path integrals involving a scattering power of body tissue. 
This work provides a simple model for such demanding applications. There is an 
approximate linearity between RMS end-point displacement and range of incident 
particles in water, empirically found in measurements and detailed calculations. 
This fact was translated into a simple linear formula, from which the scattering 
power that is only inversely proportional to residual range was derived. The 
simplicity enabled analytical formulation for ions stopping in water, which was 
designed to be equivalent with the extended Highland model and agreed with 
measurements within 2% or 0.02 cm in RMS displacement. The simplicity will also 
improve the efficiency of numerical path integrals in the presence of heterogeneity. 



PACS numbers: 11.80.La, 29.27.Eg 

1. Introduction 

The essence of radiotherapy with protons and heavier ions lies in precise control of 
incident particles that are designed to stop in tumor volume. The targeting precision 
will be inevitably deteriorated by multiple scattering in beam modifiers and patient 
body and such effects must be accurately handled for dose calculations in treatment 
planning. On the other hand, simplicity and efficiency are also essential in clinical 
practice and there has always been need for a computational method that balances 
all these demanding and conflicting requirements. 

Fermi and then Eyges (1948) developed a general theory for charged particles 
that undergo energy loss and multiple scattering in matter. A group of particles is 
approximated as a Gaussian beam growing in space with statistical variances 
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where x is the longitudinal position and y and 6 are the projected particle position 
and angle. The original Fermi-Eyges theory adopted purely Gaussian approximation 
(Rossi and Greisen 1948) with (projected) scattering power 



where E s — m e c 2 \/2tt /a as 15.0 MeV is a constant, X is the radiation length of 
the material, and z, p, and v are the charge, the momentum, and the velocity of 
the particles. The Fermi-Rossi formula ((4]) totally ignores effects of large-angle single 
scattering (Hanson et al 1951) and was found to be inaccurate (Wong et al 1990). 

Based on formulations by Highland (1975, 1979) and Gottschalk et al (1993), 
Kanematsu (2008b) proposed a scattering power with correction for the single- 
scattering effect, although within the Gaussian approximation, 

r = o.97o(i + ^V 1 + ^)S^y, (5) 
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where I = J Q X dx' / Xq(x') is the radiative path length. Although it would be difficult to 
calculate integrals HI)-© because of the embedded integral in the £ terms, Kanematsu 
(2008b) further derived an approximate formula for the RMS displacement of incident 
ions at the end point in homogeneous matter as 



* v o (Ro) = ( ~Y~ >.816 - 0.082 ln^£f-^y% /-£L (6) 

where m/m p is the ion mass in units of the proton mass, Rq is the expected in- water 
range on the incidence, ps is the stopping-power ratio of the matter relative to water, 
and k — 1.08 and A = 4.67 x 10~ 4 are constants. 

Despite the complex involvement of variable Rq in ([6J , Kanematsu (2008b) found 
the o~yQ-Ro relation for ions in water to be very linear. In fact, Preston and Kohler 
of Harvard Cyclotron Laboratory knew the linear relation and derived universal 
curve a y (x)/a y o(Ro) — \J"ix 2 R — 2xr ~ 2(1 — xr) 2 ln(l — xr) with xr — psx/Ro for 
relative growth of RMS displacement a y — (y 2 ) 1 ^ 2 in homogeneous matter in an 
unpublished work in 1968. Starting with the empirical linear relation, this work is 
aimed to develop a simple and general multiple-scattering model to improve efficiency 
of numerical heterogeneity handling and to enable further analytical beam modeling. 

2. Materials and methods 

2.1. Linear- displacement model 

Linear approximation o~ y o oc Ro for homogeneous systems greatly simplifies ([6]) to 



m \ -0-46 , y p 
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a y0 (Ro) = 0.0224 z" 08 ( — ) W^^, (7) 
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where Xq w — 36.08 cm is the radiation length of water, \J Xq w / (p^Xq) is the 
scattering/stopping ratio of the material relative to water, and Rq/ps is the 
geometrical range. Equation was calibrated to ([6]) for water at Rq = Xq- 
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2.2. Formulation of new scattering power 

Equation y 2 = a 2 at the end point x = Ro/ps associates $5$ and © as 

[ R0 *T*« = 0.0224^ ,-o." (™) -° 92 (*L V , (8) 
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to lead to another scattering power 

T = / mz ^i / m ^(l.00xl0- 3 )z-°- 16 ( — ) , (9) 
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where f mz is the particle-type-dependent factor. The scattering power is inherently 
applicable to any heterogeneous system by numerical integral of ([T]H((3]). 



2.3. Comparison of models and measurements 

We examined these Fermi- Rossi ((3]), extended Highland (O, and linear-displacement 
^ models with unpublished measurements by Phillips (Hollmark et al 2004) , those by 
Preston and Kohler (Kanematsu 2008b), and Moliere-Hanson calculations by Deasy 
(1998). We took growths of RMS displacement cr y (x) with depth x for R = 29.4 cm 
protons, i?o = 29.4 cm helium ions, and i?o = 29.7 cm carbon ions in water and RMS 
end-point displacements (J y o(Ro) for them with varied incident range Rq. 



2.4- Formulas for homogeneous systems 

For a point mono-directional ion beam with in-water range Rq incident into 
homogeneous matter with constant ps and Xq, equations dH) — dHJ) are analytically 
integrated to 

^ = / m2 ^ln^, (10) 
p s X a R 
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as a function of residual range R = Rq — psx at distance x. Equation (|12p in fact 
reduces to the universal curve by Preston and Kohler. 

A radiation field at a given x position can be effectively or virtually modeled 
with a source, ignoring the matter (ICRU-35 1984). The effective extended source 
is at x c — x — y9(x) / 9 2 (x), where y 2 would be minimum in vacuum. The virtual 
point source is at i v = x — y 2 (x)/y6(x), from which radiating particles would 
form a field of equivalent divergence. Similarly, the effective scattering point is at 
x s = x — [y 2 {x) /9 2 {x)] 1 / 2 , at which a point-like scattering would cause equivalent 
RMS angle and displacement (Gottschalk et al 1993). 



3. Results 



Figure [2(a) shows RMS-displacement growths <J y (x) and figure QJb) shows RMS 
end-point displacements <t v q(Ro). The present model was virtually identical to the 
extended Highland model with deviations from measurements or Moliere-Hanson 
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Figure 1. (a) RMS displacement cr y with depth x and (b) RMS end-point 
displacement a y o with incident range Ro for protons, helium ions, and carbon 
ions in water, calculated with the Fermi- Rossi (dashed), extended Highland 
(solid), linear-displacement (dotted— overlapping with the solid) models along with 
a Moliere-Hanson calculation by Deasy (o) and measurements by Phillips (x, o, 
+) and by Preston and Kohler (A), (c) End-point lateral shape of a proton 
pencil beam in water (Ro = 11.4 cm) measured by Preston and Kohler (A) and 
estimated by Fermi-Rossi (dashed) and the present (solid) formulations. 
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Figure 2. Analytical curves of (a) parameters 9' 2 , y8, and y 2 and (b) relative 
distances to effective extended source (e), the virtual point source (v), and effective 
scattering point (s) as a function of normalized depth x/Ro for protons in water. 



calculations within either 2% or 0.02 cm, while the Fermi-Rossi model overestimated 
the RMS displacements by nearly 10%. Figure [TJc) shows an end-point shape of a 
pencil beam measured by Preston and Kohler for Rq = 11.4 cm protons in water 
and curves estimated by Fermi-Rossi and the present formulations with additional 
(0.168 cm) 2 to (Xy for incident beam emittance in their experiment. The Gaussian 
approximation was in fact adequate in this case. 

Figure Ufa) shows analytical variances 8 2 , yO, and y 2 growing with depth x for 
protons in water. Figure [2fb) shows the relative distances from current position x 
to the effective extended source cc e , the virtual point source x v , and the effective 
scattering point x s . They approach (x — x e )/x — > 1/2, [x — x v )/x — > 2/3, and 
(x — x s )/x — > l/v3 at the x — > limit. Increase of the scattering power with depth 
moves these points relatively closer to the current position. 
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Figure 3. Scaled RMS end-point displacement ps^yO as a function of incident in- 
water range Ho and its slope pscr y o/Ro for water (solid), pure hydrogen (dashed), 
and pure calcium (dotted) calculated with the extended Highland model. 

4. Discussion 

In application of Bragg peaks, the end-point displacement is the most important, for 
which the Gaussian approximation was valid in the proton experiment by Preston 
and Kohler. In fact, the single-scattering effect is theoretically small for a thick target 
(Hanson et al 1951). Although ions suffer nuclear interactions with resultant fragments 
that generally scatter at large angles (Matsufuji et al 2005), their contributions may 
be relatively less significant at the Bragg peaks. The present formulation would be 
thus adequate for radiotherapy 

The linear-displacement model with the Fermi- Eygcs theory has brought general 
formulas for ions, including the universal curve intuitively derived by Preston and 
Kohler without explicit formulation of the scattering power. In the present model, 
the kinematic properties are encapsulated in residual range R = R — J Q psdx' that is 
always tracked in beam transport. The ion- type dependence in scattering angle and 
displacement is simply as proportional to i// mz , which leads to 50% for 4 He, 28% for 
12 C, and 24% for 16 with respect to that of protons for a given incident range. These 
numbers coincide with detailed numerical calculations by Hollmark et al (2004). 

The linearity between end-point displacement and range observed for water is the 
basis of the present model. Its validity for general body-tissue materials is not obvious. 
We here examine water and two extreme elements hydrogen (psXq 113 cm) and 
calcium (psXq ps 14.5 cm) among major elements of body tissues (ICRU-46 1992), 
using ([6]) in the extended Highland model. Figure [3] shows their Ro-(T y o relations 
with geometrical scale correction and indicates that the linearity will generally hold 
for body-tissue elements. However, the elemental linearity may not truly warrant 
the validity for systems heterogeneous in atomic compositions. The scattering power 
([5]) only depends on the residual range that is irrelevant to the multiple scattering 
accumulated in the other upstream materials, whereas the accumulation should 
influence the single-scattering effect (Kanematsu 2008b). In other words, the present 
model implicitly assumes heterogeneity in density only. 

In the current practice of treatment planning, the patient heterogeneity is 
normally modeled with variable-density water (Kanematsu et al 2003), for which the 
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present model is rigorous with further simplified scattering power 

T = f mz ^. (13) 

The simplicity will minimize computation of the integrands in path integrals ([l])-([3]) 
for demanding dose calculations (Kanematsu et al 1998, 2006, 2008a). 

5. Conclusions 

A novel multiple-scattering model has been formulated based on the fact such that the 
RMS end-point displacement is proportional to the incident range in water. The model 
was designed to be equivalent with the extended Highland model for stopping ions in 
water and agreed with measurements within 2% or 0.02 cm in RMS displacement. 

The resultant scattering-power formula that is only inversely proportional to 
residual range is much simpler than former formulations and can be used in the 
framework of the Fermi-Eyges theory for Gaussian-beam transport in tissue-like 
matter. The simplicity enables analytical beam modeling for homogeneous systems 
and improves efficiency of numerical path integrals for heterogeneous systems. The 
present model is ideal for demanding dose calculations in treatment planning of heavy- 
charged-particle radiotherapy. 
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